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ABSTRACT 

We present theoretical modelling of the electron distribution produced by annihilating 
neutralino dark matter in dwarf spheroidal galaxies (dSphs). In particular, we follow 
up the idea of Colafrancesco (2004) and find that such electrons distort the cosmic 
microwave background (CMB) by the Sunyaev-Zeldovich effect. For an assumed neu- 
tralino mass of 10 GeV and beam size of 1", the SZ temperature decrement is of the 
order of nano-Kelvin for dSph models with a soft core. By contrast, it is of the order 
of micro-Kelvin for the strongly cusped dSph models favoured by some cosmological 
simulations. Although this is out of reach of current instruments, it may well be de- 
tectable by future mm telescopes, such as ALMA. We also show that the upscattered 
CMB photons have energies within reach of upcoming X-ray observatories, but that 
the flux of such photons is too small to be detectable soon. Nonetheless, we conclude 
that searching for the dark matter induced Sunyaev-Zeldovich effect is a promising 
way of constraining the dark distribution in dSphs, especially if the particles are light. 

Key words: galaxies: dwarf - intergalactic medium ~ cosmic microwave background 
- X-rays: galaxies - cosmology:dark matter 



1 INTRODUCTION 

Dwarf spheroidal galaxies (dSphs) are important probes of 
dark matter. They are among the highest mass-to-light sys- 
tems known, and the dynamics of their sparse stellar pop- 
ulations is governed by the dominant dark matter distri- 
bution. In addition, no emission has been detected from 
dSphs in wavebands other t han the optical, indicating a 
lack of internal dus t or gas jFomalont fc Geldzahle^ll979l: 
iBonanos et"alll2004f) . 

Here, we consider the distortion of the cosmic mi- 
crowave background (CMB) by the non-thermal popula- 
tion of secondary electrons generated by dark matter an- 
nihilation. This is an example of the Sunyaev-Zeldovich 
(SZ) effect iZeldovich fc Sunvaevl ll969t) . Since we are deal- 
ing with electrons produced from dark m atter annihilation, 
we wr ite t he distortion as the dSZ ef fect. lEnsslin fc Kaisen 
i2000l) and lColafrancesco et al.l i2003l) calculated the signal 
expected from the SZ effect of a relativistic plasma, while 
fColafra ncesco. 120041 determined the dSZ effect in galaxy 
clusters. 



* E-mail: tlc26@mrao.cam.ac.uk, 
cola@mporzio.astro.it 



nwe@ast.cam.ac.uk. 



As first pointed out in IColafrancescol i2004l) . dSphs are 
attractive targets because they have very high mass-to-light 
ratios and because they have few contaminants. In partic- 
ular, they have little or no internal magnetic field, so it 
is not possible for synchrotron emission (from annhilation 
electrons or otherwise) to contaminate the dSZ effect. Since 
dSphs are believed to be almost devoid of interstellar gas, 
other mechanisms such as HI or CO line emission will not 
be important. CMB distortions are therefore a rather clean 
method to detect the annihilation signature. We therefore 
use existing models of the dark matter distribution in dSphs, 
and a possible form of the energy spectrum of electrons pro- 
duced by dark matter annihilation, to derive the tempera- 
ture change in the CMB. We focus on dark matter models 
explicitly constrained by observations, rather than simula- 
tions. 

For our predictions, we assume that the cold dark 
matter particle is the lightest suporsymmctric particle, the 
neutralino (Ju ngman. K amiorikowski & Gricst 1996). Cur- 
rent limits on the neutralino mass and centre-of-mass 
velocity-av eraged cross-section (aV ) A hav e been reviewed 
recently bv lBertone. Hooper fc Sili3 l|200^. Based on these 
results, and considerations of current or upcoming exper- 
iments, we investigate the SUSY parameters {a-V)A = 
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Table 1. List of CHM model parameters for the Draco dSph. Val- 
ues for rt in parentheses are for NFW models of the Milky Way, 
as opposed to isothermal power-law models which are without 
brackets. 



7 A X rs/kpc n/kpc A/(rt)/lO**M0 



0.5 


2.3 


0.32 


6.6 (1.5) 


5.5 


1.0 


3.3 


0.62 


7.0 (1.6) 


6.6 


1.5 


2.9 


1.0 


6.5 (1.5) 


5.5 



Table 2. List of CPL model parameters for the Draco dSph. Val- 
ues for rt in parentheses are for NFW models of the Milky Way, 
as opposed to isothermal power-law models which are without 
brackets. 



a 


■Ua/km s 


rc/kpc 


n/kpc 


M{rt)/IQ^MQ 


0.2 


24.7 


0.25 


6.2 (1.3) 


4.6 





22.9 


0.23 


7.8 (1.4) 


9.5 


-0.2 


20.9 


0.21 


10.1 (1.6) 


22.43 



10~^^ cm^s~^ and = 10 GeV. Such low mass particles 
may provide a sizeable contribution to the matte r density in 
the Universe llBottino. Fornengo fc ScopeJl20o'3) . and hence 
are worthy of consideration. However, some dark matter can- 
didates - such as the neutralino in the most commonly stud- 
ied minimal supersymmetric models or the lightest Kaluza- 
Klein particle - must be more massive than this. The value 
assumed for {aV) a is consistent with the expected relic den- 
sity in a Universe with = 0.3 and Hq — 70kms~^Mpc~^ 
JColafrancesco fc Meldliool .We derive the dependency of 
dSZ signal on {aV) a and M^, and show that the bright- 
ness temperature decrement ATb oc {aV) aM~''^ . We ini- 
tially calculate the expected signal for this optimistic choice 
of dark matter parameters. 

The format of the paper is as follows. In Section |5| we 
discuss our dSph dark halo models, based on current work 
in the literature. We then discuss possible products of dark 
matter particle annihilation in Section|3 Observational con- 
sequences of such annihilation events are presented in Sec- 
tion 2] We conclude in Section |S] 



2 DWARF SPHEROIDAL MODELS 

With spheric al symmetry of the dark halo a ssumed, we use 
the results of lEvans. Ferrer fc SarkaJ i2004) . who fit obser- 
vational data o n the Draco dSph (cur rently orbiting the 
Milky Way) from lWilkinson et al. | l|2004'l to two sets of mod- 
els via the Jeans eauatioiT lBinng^g Treinain& .1987il . 



2.1 Cusped Halo Models 

Cusped halo models (CH Ms)are favoured by numerical sim- 
ulations, as reported b y iNava rro. Frenk fc Whit^ (1199711 
(NFW) and iMoore et al. (1998)). Using arguments based 
on the survi vability of kinematically cold substructure, 
iKlevna et all (.20011 argued against a cusped halo for at least 
the Ursa Minor dSph. For the rest of the dSphs, cusped halos 
remain viable. We therefore consider number density radial 
profiles of the form 



n (f) = noa (f) , 



(1) 



where no = Ar~ /M^, and a (f ) = f^^ {I + fy~ , with 
r = r/vs- Table gives the cusp slope 7, the scale radius 
Ts, the tidal radius rt, and the overall normalisation A. The 
model is truncated at rt, whose value depends on the Milky 
Way dark halo model. 

To address the problem of the divergent cen- 
tral density in this mode l, w e use the a rguments of 
iBlasi fc Colafrancescd ^1999^ and lTvleJ J2002ll . The density 



profile is truncated at a radius r,nin, where we assume the 
dark matter annihilation rate matches the collapse timescale 
of the cusp. With this assumption, a small constant density 
core is created with radius rmin 



rt{aV)T ( 



\ AtvGM^ 



1/4 



(2) 



where nmin is the number density at the location of the 
constant density region. At small radii, Eq.(0 reduces to 



n(r) ~ AM- 



^, and hence we find: 



^1+7/4 



rt{aVy/ 



1/2 



7-3 N 1/4 



(AiLLX 



(3) 



Typically, rmin ~ lO^'^m, and for such a small value, it is 
unlikely that tidal forces (from M31 or the Milky Way) 
could disrupt the central cusp. It should be noted that the 
Moore profile (7 = 1.5) represents an extreme case for the 
inner slope of the cusp. Recent numerical simulation s (e.g. 
iNavarro et al l2004l:lDiemand. Moore fc Stadel 120041) point 
towards a milder density slope of 7 ~ 1.1. However, we re- 
tain the Moore profile here as the models satisfy the observa- 
tional constraints from stellar radial velocities. In addition, 
this profile allows an upper limit of the magnitude of the 
dSZ effect, enabling a broad yet well motivated parameter 
space to be investigated. 



2.2 Cored Power-Law Models 

The second family of dark halo profiles s tudied are the cored 
power- law (CPL) models jEvanslll99J) . Again, these sat- 
isfy the Draco velocity dispersion observations, and take the 
form 



n (f) = nob (f ) , 



(4) 



where no = v^/'i-jiGr'^. The radial dependence is given by 



6(f) 



3-f f2(l-a) 



(5) 



(l+f2)2+«/2' 

where r = r/rc. Table |5|gives the slope a, the core radius r^, 
the tidal radius vt and the velocity scale Va- The dSph dark 
halo extends out to rt, the value of which again depends 
on the model (isothermal power-law or NFW) used for the 
Milky Way. 



3 DARK MATTER ANNIHILATION 

Having established a set of representative dark matter halo 
models for dSphs, which embrace a range of possible struc- 
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tures, we determine the decay products of annihilating neu- 
tralinos. Electrons produced in this fashion have an as- 
sociated cooling function b (E) , comprised of synchrotron 
losses, inverse Comption scattering (ICS) and Coulomb 
losses JColafrancescoll2004l) : 



HE) 



dE 
'dt 



dE 
'dt 



+ 



f dE 
Vdl 



GeV s" 



(6) 



At el ectron energies E > ISOMeV ijBlasi fc Colafrancescol 
Il999 h. this function is dominated by the first two terms, so 
the Coulomb term may safely be dropped from the cool- 
ing function. In addition, there is no evidence in favour of 
signifi cant magnetic fields in dSphs iFomalont fc GeldzahleJ 
I1979I) . We therefore also drop the sy nchrotron cooling term, 
leaving l|Colafrancesco fc Melel|200lfl 



b{E) = 2.5 X 10" 



'^"^ xE^ GeV s" 



(7) 



Energy losses are efficient in the dSph halo, and dark matter 
annihilation from infalling material continuously refills the 
electron spectrum. In a first approximation, it seems reason- 
able to neglect difi'usive effects, as dark matter annihilation 
replenishes the electron spectrum efficiently, at least in the 
central parts of the dSph. 

Furthermore, we consider the dark matter to maintain 
a constant radial profile on the timescale of dark matter an- 
nihilation. This means that any time-variation of the elec- 
tron distribution is negligible. The source function of elec- 
trons, Qe {E, r), generated in this fashion therefore obeys the 
stationary diffusion-loss equation (see, for example, Longair 
1994) 



_d_ 

dE 



dUe {E, r) 
dEe 



HE) 



-Qe {E, 



(8) 



which we integrate to obtain the annihilation-produced elec- 
tron energy distribution dUe {E,r) /dEe'. 



1 



HE) 



Qe{E,r)dE. 



(9) 



The source function itself results from annihilation 
products of neu tralino colli sions. In this calculation, we fol- 
low the work of lTvlerl fcoo'sl) . who used the HiU (1983) for- 
mula to determine the densities of particles produced by 
dark matter annihilation. Quark pairs and their subsequent 
fragmentation lead to pions as the main annihilation prod- 
ucts. The primary decay particles are neutral pions, which 
decay to gamma rays, and charged pions, which decay as 



and TV 



M I'm- 



The muons then decay to electrons via 

/i^ e'^Ufj^Ue and fi^ e 



(10) 



(11) 



The number spectrum of electrons from a single xx annihi- 
lation is then given by: 



dEe 



where r = (m^/m^)^, and the charged pion multiplicity per 
annihilation event is 



4 15 



\ -3/2 / E \^ 

3 16M^ V j V'lt) ' ^^^^ 



where the factor of 4/3 accounts for the fact that annihila- 
tion electrons are only produced by charged pions, and that 
quarks (which eventually decay to charged pions) are pro- 
duced in pairs. The number spectrum of muons produced 
per charged pion decay is 



and 



dm 

dEe 



2 

e: 



dE^ 



E„ ml 



2 

+ 3 



(14) 



(15) 



is the number spectrum of electrons per muon decay. After 
some algebra, with the above forms for the decay product 
energy spectrum, Ea. l|12|l has an analytic solution: 



dNe 



15 



ClZ 



-3/2 



+ 



-1/2 , I 1/2 . 2 , 3 

+ C3 -f C4,Z -f- C52: + CaZ 



(16) 



in units of GeV ^, where z = Ee/M-^. The coefficients 
Ci are (0.1039,-1.2218,2.4800,-1.5406,0.2205,-0.04197). 
Scaling this expression, we arrive at the source function 



( 



3 8 m| - m2 V GeV 



10- 



-3/2 , -1/2 

ClZ ' + C2Z ' 



I 1/2 I 2 I 3 

-I- ca + C4Z ' + C5Z + C(iZ 



(17) 



where the appropriate units are GeV~^cm~ s~^. This needs 
to be multiplied by a further factor of 2 to account for the 
contributions of electrons and positrons. Using this expres- 
sion, and Eq.@, we arrive at the equilibrium electron energy 
distribution 



drie 
dK 



■ 4.6848 X 10" 



10-26cm^ 



y 

GeV J 



diz 



5/2 



+ d2Z 



-3/2_ 



d'iZ ^ + d. 



,4Z 



-1/2 



-f dsz + dez 



(18) 



in units of GeV~ cm~ . In this expres- 

sion, the radial dependence of due/dEe arises 
in the factors. The new coefficients di — 

(0.2078, 2.4436, -2.4800, 1.0271, -0.07333, 0.01049). Equa- 
tions 1171 1 and II8I I are displayed in Figs0and|21respectively. 

For completeness, in Fig we also display two other 
possib le source functions presented in lColafrancesco fc Meld 
(l200lfl . These source functions arise from fermion- 
dominated annihilation in the first instance (denoted CM 
XX ~^ ff)i and gauge boson dominated annihilation (CM 
XX ~* WW) in the second. At the crucial low energy end 
of the spectrum, the three separate source functions ex- 
hibit amplitudes within 1.5 orders of magnitude, and simi- 
lar power-law slopes. This implies that our calculations may 
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Hill 

CM xx-ff 
CM XX^WW - 




Figure 1. Electron source function, as described by Eg 1171 (Ae- 
noted Hill), assuming = 100 GeV, (aV)A = IQ-^*^ cm^s-l 
and rijj- = lcm~^ . Also shown are the source functions described 
in Appendix A o l lColafrancesco fc Mel j <200ll) , i.e. electrons pro- 
duced in fermion-dominated annihilation XX ~* //i ^^'^ gS'Uge- 
boson dominated annihilation, xX ~* WW. 



> 



— ^ — I.J — I — VJ.J u 

logio(z) 

Figure 2. Equilibrium electron energy distribution, calculated 
from Eq[I21 using AI^ = 10 GeV, {ctV)a = 10" 
?iv = lcm~ 



-3 



for Qe from gener al electron sources are possib le, such as a 
double power-law JColafrancesco fc Meldliool . these again 
have negative slopes and therefore it is the lowest electron 
energies that are most significant. 

From Eq l|TSl, we may calculate the number density of 
electrons in the halo: 

" dne 



■dE, 



The limits are chosen to be O.OIM^ < E^ < M 



(19) 



based on 



the analysis of iKamionkowksi fc Turner! jigOlT) . who calcu- 
lated the positron source function in a similar manner to 
that described above. We will see that the upper limit is 
somewhat irrelevant, as the source function falls rapidly to 
zero as the electron energy approaches the neutralino rest 
mass energy. 

Evaluating this leads to 



ne = 8.09 X 10^ 



10-26cm» s-i J V GeV 



(20) 

in units of cm ^. 

We now represent Eq 1)18^ in terms of its momentum 
spectrum, as this quantity will be used later to calculate the 
frequency spectrum of upscattered CMB photons. Specifi- 
cally, we write 



^ = (r) fe (p) , 



(21) 



where p — /3e7e is the normalised electron momentum, and 
fe (p) has the property 



/ dp/, (p) = 1. 

Jo 



(22) 



In the case of dark matter annihilation, f3e — 1.0 to a very 
good approximation. Furthermore, 7e = Ee/meC^, which 
will be in the range 19. GM^ < 7e < 1960Mj(., corresponding 
to energies in the range 0.01 — 1.0 M;^/GeV. Eq II21II can 
be written explicitly by retaining the power law terms in 
Eq. 1181 and their associated 'weights' di, and introducing a 
normalising factor A = 1/172.77. The momentum spectrum 
is then written 

/e(p)=A^ d,f^V (23) 



where 
have 



(2.5,1.5,1.0,0.5,-1.0,-2.0). Finally then, we 



= 1,6 



(24) 



increase or decrease by approximately one order of magni- 
tude in either direction, depending o n the precise details of 
the so urce function. As the results of lColafrancesco fc Melel 
l|200l|) wer e derived in an independent manner to those of 
iHill I il983l) , our results may also hold for neutralino compo- 
sitions other than those considered here. The crucial quan- 
tity for the magnitude of the dSZ effect is the lower limit of 
electron energy, which can increase rapidly for any power- 
law source function Qe ~ where /3 > is a generic 
slope parameter. It is highly unlikely on energetic grounds 
that electrons produced by annihilating dark matter can 
have Qe rising with z. Although more complicated forms 



4 OBSERVABLES 

We now calculate the magnitude of dSZ radio emission, and 
the flux of up-scattered CMB photons, using the model pa- 
rameters for Draco as described in Section |2| The compu- 
tations use an approximate method based on the work of 
Ensslin & Kaiser ( 2000) . This holds to first order in the elec- 
tron optical depth, neglects multiple scatterings and is vahd 
for a single electron population only. Throughout, we em- 
ploy units where {aV)A is measured in 10~^^cm^s~^, and 
in GeV. Distances such as and rc are taken in kpc, 
and nx,o is in cm~^. 
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4.1 The Sunyaev-Zeldovich EfFect 

The SZ effect is caused by inverse Compton scattering of 
CMB photons off energetic electrons. On average, the pho- 
tons gain energy, shifting their spectrum to higher fre- 
quencies, and causing a distortion in the CMB radiation 
field. This is commonly characterised by the Compton y- 
parameter, the line-of-sight integral of gas pressure through 
the electron cloud: 



-Al. 



(25) 



We may also compute the integrated y-parameter, Y , over 
the solid angle Q. of the dSph: 



Y = 



(26) 



An integral over solid angle can be written as df2 = AA/r\, 
hence we can equivalently write 



Y = 



PAV, 



(27) 



where rn is the heliocentric distance to the dSph (assumed to 
be 80kpc for Draco. For the ultra-relativistic gas considered 
here, where Ee ~ GeV, we apply the relationship between 
pressure and energy density Pe = e/3. The (domi nant) ki- 
netic energy density e is obtained from (lEnsslin fc Kaiser! 
I2OOOI) 

€ = n,{r) dp/4p)(vTTp2-l)meC^ (28) 

where again the upper limit is not too important, as the 
source function falls rapidly to zero as the electron energy 
approaches the neutralino rest mass energy. We already 
showed that p > 19.5695Mx, and so even for a low neu- 
tralino mass of lOGeV, / » 1. In this case, we can simplify 
Eq. to 



= ne{r)meC / dpfe{p)p. 

Jo.OlM^/m^ 



(29) 



The pressure is then evaluated using Eqs ^ and fM . 
leading to: 



Pe = 1.585 X 10"^^n^((jy>A J m-^ 



(30) 



Note that the neutralino mass is present only in the number 
density term here. 

We may now calculate the integrated Y parameter, us- 
ing this result and Eq. I27II . for each dSph model in Section 

m 



Y = 3.974 X 10-^rlrJ^^nl^o{(^V)A 



i'^a' (f) df, (31) 



where Ts and a (f) should be replaced with and h (f) for 
the CPL models. 

To derive temperature shifts, we will work in terms of 
the mean Compton parameter averaged over the dSph, i.e. 



Y 



(32) 



where f2 is the angular extent over which the Y parameter is 
averaged in Eq. 12611 . Once converted to temperature units, 
y measures the temperature decrement inside a telescope 
beam of angular size f2. We choose Q. to take three values, 



first to match the angular size of the whole dSph, secondly 
within a 1' beam, and finally within a 1" beam. 

For the assumed neutralino mass, annihilation electrons 
are always ultra-relativistic. Since we expect these particles 
to have energies ~ GeV, the effect of such a non-thermal 
electron population is to completely remove photons from 
the spectral range of the CMB. The problem is one of elec- 
tron number density: a low neutralino mass increases the 
number density of electrons as M^, which raises the scatter- 
ing probability accordingly. It is clear that the relativistic SZ 
signature really measures the electron number in the dSph. 
Since there is a one-to-one correspondence between the num- 
ber of electrons and the number of annihilating neutralinos, 
y is a direct measure of the dSph mass. 



J^.l.l The Intensity Shift 

The fractional change in the CMB intensity field is 



Si (x) = 



A/ {x) 



(33) 



where x — hu/kbTcMB is the dimensionless frequency, and 
lo ~ 2 (/cTcmb)^ / (hc)^ . Contributions to such a distortion 
are written as a product of a spectral function, g (x), and the 
Compton parameter y. The product g (x) y usually refers to 
the thermal electron population in clusters of galaxies; for 
the relativistic gas considered here, we write g (x) to make 
the distinction explicit. The fractional distortion averaged 
over the whole dSph is therefore written as (iRaphaelL. 1 99 sl) 



Si (x) =g(x)y. 



(34) 



There are two contributions to g{x). First, photons are re- 
moved from the infinitesimal frequency band x,x ~\- /S.x by 
collisions with the ultra-relativistic electrons. This contri- 
bution is written —i [x). The second effect is the photons 
scattered into this band from lower frequencies, which is 
written j {x). Scattered CMB photons have their frequency 
increased (on average) by a factor of 47^/3 — 1/3, so CMB 
photons are up-scattered to the X-ray regime. Therefore, we 
are interested in two distinct frequency bands - that near 
X ~ 2.5 in the radio corresponding to i{x), and the band 
close to a; ~ lO" in X-rays, described by j {x). 

Essentially no photons are scattered into the radio fre- 
quency band under consideration from lower frequencies. 
Photons are simply removed from the spectrum, causing a 
decrease in the number of photons, and thus a corresponding 
decrease in the specific intensity. The factor i (x) therefore 
has the spectral form of the CMB i (x) = x^ / (e^ — 1), which 
is maximal at x = 2.82. 

The spectral fact or 
by jEnssUn fc Kaiseill2000t) 



5 (a;) 



thus 



given 



g (x) = [j (x) - i (x)] 



kbT, 



(35) 



We define the pseudo-temperature ktTe, as the ratio of the 
gas pressure P^ to the electron number density n^. For a 
thermal electron distribution, this is equal to the thermody- 
namic temperature. In this case, we have 



kbTe = — = 0.0122426 



GeV 



GeV. 



(36) 
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Table 3. List of y parameters, fractional intensity shifts Si, and temperature decrements AT for both sets of dSph models, assuming 
= 10 GeV and {(tV)a = lO"^** cm^ s'^ All quantities are averaged over the angular extent of the dSph. Bracketed values correspond 
to the different Milky Way models as described in Table [Tl A heliocentric distance to Draco of 80kpc is assumed. For the intensity and 
temperature shifts, we use a frequency of a; = 0.616 i.e. 35GHz. 



7 


5/10-" 




AT/lO-i^K 


0.5 


7.019 (1.338x10^) 


8.026 (1.530x10^) 


-5.973 (-1.103x10^) 


1.0 


7.095 (1.336x10^) 


8.114 (1.528x10^) 


-6.038 (-1.137x10^) 


1.5 


4.514x10^ (9.365x10^) 


5.162x10^ (1.071x10^) 


-3.841x10^ (-7.970x10^) 


a 




5i/10-i3 


AT/10-13/K 


0.2 


4.045 (9.040x10^) 


4.626 (1.034x10^) 


-3.442 (-7.693x101) 


0.0 


2.970 (8.599x10^) 


3.396 (9.833x101) 


-2.527 (-7.318x101) 


-0.2 


2.250 (7.586x101) 


2.573 (8.676x10^) 


-1.915 (-6.456x101) 



Using Eqns l|84|l and H35|l . the fractional distortion in the 
CMB at radio wavelengths, integrated over the dSph, is 
therefore: 



Si {x) = ~i {x) ^^^y. 



(37) 



The approximation of drop ping j{x) is valied provided x < 



ine approximation ot drop 
10 jEnsslin fc Kaiseill2000l) . 



4.1.2 The Temperature Shift 

The above result can equally be expressed as a temperature 
shift. Expressing the result in this manner has the advan- 
tage of a direct comparison to typical CMB telescope noise 
temperatures. 



AT : 



I dl (x) 
I dT 



AI 



(38) 



gives the expected temperature shift in the CMB, where 
I (x) = loi (x). The partial derivative is 



dl{x) 



(39) 



dT " (e- - 1)2 ' 
and in conjunction with Eqs. 1331 . 1371 and 13811 . we have 



AT (a;) 



1 m„c 



■yTcMB. 



(40) 



Table shows the results for the fiducial case of = 
10 GeV and {aV)A = 10"^® cm'^ s~i, assuming a helio- 
centric distance ru = 80 kpc for Draco and with a beam size 
that matches the angular size of the dSph. We assume an 
observing frequency of 35GHz [x — 0.616, close to Band 1 of 
the forthcoming Atacama Large Millimetre Array (ALMA)). 
While such y values are prohibitively low for current or near- 
future exp eriments - for exa mple, the upcoming South Pole 
Telescope ijRuhl et al.ll2004^ will only achieve noise temper- 
atures of ~ /iK on arcminute scales. However, noise tem- 
peratures of order ~ 10~^ K may well be reached by a fu- 
ture generation of radio/sub-mm telescopes, such as ALMA. 
Finally, Tables EE] show the same quantities, but for beam 
sizes of 1', 1" and 0.1", encompassing a range of target spec- 
ifications for ALMA. 

It is instructive to compare the results presented here 
with those of Colafrancesco (2004) on clusters of galaxies. 



Although clusters are considerably larger objects, their cen- 
tral dark matter number density is at least a factor of 10 
lower than in the dSph models considered here. Since the 
electron density is proportional to n^, we therefore have a 
relative increase in dSZ pressure of at least a factor of 100 in 
dSphs. Using Eq l|25^ . it is also clear that the integral along 
the line of sight will be proportional to the scale length of 
the object, of order 10^ kpc for a cluster and lO^i kpc for 
dSphs. Hence, these two competing factors almost cancel, so 
the surface brightness for the dSZ effect is roughly equal for 
clusters and dSphs. 

In the most optimistic cases, the dSZ effect is within 
the grasp of ALMA. For example, using the 7 = 1.5 cusped 
model at lOOGHz (i.e. ALMA Band 3) and converting to 
brightness temperatures Tb rather than the thermodynamic 
temperatures quoted in the Tables, the ALMA sensitivity 
calculator for 64 dishes with resolution 1" gives a signal- 
to-noise of unity in two hours. We present theoretical con- 
straints in the Af^ — Tb plane in Fig. |21 The plot displays 
the brightness temperature for cusped models observed in a 
1" beam with 24hr and 1 yr integration times at lOOGHz. 
This plane contains a substantial portion of the parameter 
space considered viable in particle physics. It is only for the 
most extreme cusp (7 = 1.5) and low neutralino masses that 
a detectable signal is predicted. 

However, we caution that the situation examined here 
may be over-simplified for the cusped models. If such den- 
sity profiles really are followed to small radii, the Coulomb 
term neglected in Eq. Qmay become important. This effect 
requires a more thorough treatment to assess its impact - 
although extra cooling would lower the average electron en- 
ergy, it would also drive the electrons closer to a thermal 
equilibrium. The spectral factor would in turn be closer to 
the thermal form, which tends to have larger values at the 
frequencies considered here. The CMB temperature shifts 
presented here may therefore be an underestimate at the 
cusp centre. 

In addition, it is possible that the electron distribution 
may deviate from the dark matter cusp structure once the 
Coulomb term becomes dominant. In that case, the electron 
density distribution may become less steep, and more like 
a cored profile. On the one hand, this means that the peak 
signal will be lower. On the other, the signal averaged over 
a larger beam increases, as the effective cusp more closely 
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ALMA 1 year 



7=1.5 



7=1.0 



7=0.5 



10 100 1000 

M^/GeV 

Figure 3. — Tb plane for three cusped dSph profiles, with an 
assumed NFW profile for the Milky Way (isothermal power-law 
models give almost identical results). We assume a 1" beam at 
lOOGHz, and display the brightness temperature sensitivity for 
24 hr and 1 yr ALMA observation. 



matches the beam size. The tabulated temperature decre- 
ments increase with decreasing resolution precisely for this 
reason, and so a 'smoother' cusp is a better match to the 
angular scales of the telescopes mentioned above. 



4.2 Upscattered Photons 

We now consider explicitly the spectrum and flux of upscat- 
tered CMB photons, adopting the flrst-ordcr, approximate 
approach described iu iEnsslin fc Kaiser. (.2000.') . 



4-2.1 Frequency Speetrum 

The scattered spectrum of photons, j [x), can be expressed 
as 



j {x) = / dtp {t) i {x/t) , 

JO 



(41) 



where the photon redistribution function P (t) gives the 
probability of a photon being scattered to a frequency t 
times greater than its original frequency. If the electron 
momentum spectrum is fe (p) dp, the photon redistribution 
function is 



P{t)= / dp,U(p)P(t-p). 



(42) 



In this expression, P{t;p) is the redistribution function for 
a monoenergetic electron distribution. This has an analytic 
form 



32pH 

3{l + t) [- 3-H3p^ +/ 
3 + 2p2 



2p 



(2arcsinh (p) — | In (t) 



(43) 



with the condition that P {t;a,pi,p2) = if | In (t) | > 
2arcsinh (p2). We may apply this formalism to the nor- 
malised electron momentum spectrum in Eq. II23II . For such 
high electron energies, it is more convenient to express 
Eq 1431 in terms of the logarithmic frequency shift s = In (t), 
in which case: 



P (s; p) = P (e*; p) e^ds. 



(44) 



P{t) is plotted in Figure^ for A'4 = lOGeV and {aV)A = 
10~'^^cm'' s~^ as before. Having computed this function, we 
calculate j (x) via Eq. 14111 . In Fig|3 we plot \j (x) — i {x) |, 
which clearly displays the prominent emission in the radio 
and X-ray bands. The two regions are distinctly separated, 
with i (x) constrained to radio frequencies, while j (x) dom- 
inating in X-rays. There is virtually no overlap between the 
two. In addition, j (x) peaks close to a: ~ 2.0 x 10^, corre- 
sponding to p ~ 500, as might be expected on the basis of 
the power- law nature of Eqn. 12311 . and the expected average 
electron momentum p = 47^/3 — 1/3. 

4.2.2 X-ray Flux 

The X-ray intensity produced by the upscattered photons 
can be calculated from Eq. 1341 . In this instance, we replace 
g (x) j (x), and then the X-ray flux density is 



Jxray (x) = loj (x) Y 



(45) 



where we integrate over the solid angle of the whole dSph. 

We consider the X-ray emission integrated over a uni- 
form efficiency energy band 0.1 — lOkeV (corresponding to 
xi = 4.247 X 10^ and X2 = 4.247 x lO'^), as an approximation 
to the current X-ray satellites Chandra and XMM-Newton. 
Dividing Eq. 14511 by x and integrating over the bandpass 
yields the X-ray photon flux 



-Ftot {xi,X2) 



IqY 

kTcMB 



dx. 



(46) 



For Mx = 10 GeV and {c7V)a = 10"^*^ cm^ s"\ the integral 



above evaluates to 



dx = 0.414. 



(47) 



The results for each of the models considered previously 
are listed in Table |7| The X-ray fluxes are all of order 
10~^^cm~^s~^, above current estimates of the gamma-ray 
flux from direct annihilation channels, but well below what 
is currently possible to observe with X-ray satellites. 



5 CONCLUSIONS 

In the above analysis, we have shown the following: 

• The Sunyaev-Zeldovich effect caused by secondary elec- 
trons produced from dark matter annihilation in dwarf 
galaxies (the dSZ effect) proposed bv IColafrancescJ ll2004ll 
could be measurable. The Comptonisation parameters av- 
eraged over the angular size of the dSph are y ~ lO"'^^ for 
low neutralino masses of lOGeV and {aV)A = 10~'^®cm'^s^^. 
The temperature decrement for an assumed beam size of 1" 
is of the order of milli-Kelvin for extremely cusped dSph 
halo models and a few tenths of a nano-Kelvin for cored 
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Table 4. List of y parameters for 1' beam, fractional intensity shifts Si, and temperature decrements AT for both sets of dSph models. 
The remaining dSph parameters are as recorded in Table Isl We assume = lOGeV and {aV)A = 10~^'^cm^s~^ . For the intensity and 
temperature shifts, we use a frequency of a: = 0.616 i.e. 35GHz. 



'Y 
1 


i7/10~'^ 


(5i/10""' 




AT/10~^K 


0.5 


3.209 (3.208) 


3.670 (3.590) 




-2.731 (-2.671) 


1.0 


9.619 (9.582) 


1.100x10^ (1.096X 


IQi) 


-8.186 (-8.155) 


1.5 


5.119x10^ (5.499x10^) 


-5.854x10^ (6.289X 


10^) 


-4.356x10^ (-4.680x10^) 


a 


g/10-8 


5i/10-"' 




AT/10- "VK 


0.2 


4.399 (4.397) 


5.031 (5.029) 




-3.744 (-3.742) 


0.0 


4.473 (4.470) 


5.115 (5.112) 




-3.806 (-3.804) 


-0.2 


4.418 (4.414) 


5.053 (5.048) 




-3.760 (-3.757) 



Table 5. The same as Table HI but for a 1" beam. 



7 


5/10-6 


(5j/10-3 


AT/IO-^K 


0.5 


1.066 (1.057) 


1.219 (1.209) 


-9.075 (-8.999) 


1.0 


8.102x10^ (7.251x10^) 


9.266x10^ (8.292x10^) 


-6.895x10^ (-6.171x10^) 


1.5 


1.356x10^ (1.466x10^) 


1.551x10^ (1.677x10^) 


-1.154x10^ (-1.248x10®) 


a 


g/10-8 


<5i/10-" 


AT/lO-iO/K 


0.2 


4.528 (4.471) 


5.178 (5.113) 


-3.853 (-3.805) 


0.0 


4.624 (4.548) 


5.288 (5.201) 


-3.935 (-3.871) 


-0.2 


4.599 (4.497) 


5.259 (5.143) 


-3.914 (-3.827) 



Table 7. List of integrated X-ray fluxes for Draco, in the energy 
band 0.1-lOkeV, for = 10 GeV and {aV)A = IQ-^® g-l 



7 


FxoT/lO-^^cm-^s-i 


0.5 


2.604 (2.564) 


1.0 


2.962 (2.914) 


1.5 


1.625x10^ (1.795x10^) 


a 


FTOT/10-"cm-2s-i 


0.2 


1.325 (1.301) 


0.0 


1.539 (1.436) 


-0.2 


1.955 (1.654) 



models. This may provide a definitive test between these 
competing hypotheses, if the signal for cusped model is large 
enough to be detectable by future radio telescopes. This re- 
sult holds before the noisy effects of primordial CMB, radio 
point sources, and SZ from clusters, has been taken into ac- 
count. This, however, is only a concern for cored models, 
in which the signal comes from the bulk of the dSph. Even 
then, dSphs are clean and uncontaminated objects, devoid of 
magnetic fields and gas, and mostly free from point sources 
such as supernova remnants. For cusped models, most of 
the signal comes from the very centre and so contaminating 
point sources are not a worry. 

• Our main aim here has been to demonstrate the fea- 
sibility of measuring the dSZ effect. We caution that our 
calculations make use of an approximate, single-scattering 



5 10 15 20 

s 

Figure 4. Photon redistribution function, for the momentum 
spectrum Eq. I24i . assuming M-^ = 10 GeV. 



formalism that holds good for low electron optical depth. 
We may therefore have underestimated the size of the effect 
in the very innermost regions of cusped dSph models. Fur- 
ther numerical treatments, for example using the methods 
of Colafrancesco et al. (2003), in the vicinity of dark matter 
spikes are desirable. 

• Upscattered CMB photons lie in the X-ray band, 
with the emission peak near x = 2.5 x 10^ for the neu- 
tralino mass considered here. Their integrated fluxes are 
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Table 6. The same as Table 1^ but for a 0.1" beam. 



7 


J//10-6 


(5j/10-9 


Ar/io-«K 


0.5 


1.636 (1.515) 


1.871 (1.732) 


-1.392 (-1.289) 


1.0 


1.544x103 (9.999x10^) 


1.766x10^ (1.143x10^) 


-1.314x10^ (-8.505x102) 


1.5 


7.893x10'^ (9.879x10'^) 


9.026x10" (1.130x10^) 


-6.717x10'^ (-8. 407x10") 


a 




<5j/10-io 


AT/10-i°/K 


0.2 


5.200 (4.605) 


5.946 (5.266) 


-4.249 (-3.919) 


0.0 


5.440 (4.694) 


6.221 (5.368) 


-4.630 (-3.995) 


-0.2 


5.601 (4.665) 


6.405 (5.335) 


-4.767 (-3.970) 




10 



Figure 5. logjg \j {x)—i (x) \ against logjQ (x), for the momentum 
spectrum Eq. 1241 . Note the two prominent bands, in the radio 
near x ~ 2.5 and in X-rays near x = 2.5 X 10^. 



~ 10~^^ cm~2 s~^, comparable in size to that from gamma- 
rays produced by direct annihilation channels. However, 
even next generation X-ray satellites such as Constellation- 
X, with collecting areas of ~ lOOOcm"^ will struggle to de- 
tect such a signal. 

• Our assessment of the importance of the dSZ effect is 
quite optimistic. If dark haloes are strongly cusped, then we 
conclude that the dSZ effect may be measurable in the near- 
future by telescopes like ALMA. However, if dark haloes 
are only weakly cusped, or if the dark matter particles 
are heavy (> 50 GeV), then even the most generous in- 
tegration times with ALMA may not yield a positive de- 
tection. Nonetheless, it is worth bearing in mind that here 
are some circumstances in which a larger effect may be pro- 
duced. First, the recently-discovered very dark dSph Ursa 
Major iWillman et al.ll2005l : iKlevna et alJl2005^ may be the 
first representative of the missing dark satel lites predicted 
by numerical simulations iMoore et al.|[l993) . In this case, 
there may be undetected, very dark dSphs much closer to 
us than Draco, which is beneficial as the flux received ob- 
viously varies like the inverse square of distance. Second, 
our calculations apply only to the case of the neutralino 
dark matter candidate. There are other possibi lities, in- 
cluding light (1-5 MeV) scalar dark matter (iBoehm et al. I 



120041: iHooper et alj '2004*) and Kaluza-Klein dark matter 
iBerton^^foope^;: Silk.2005k) . whose induced dSZ signals 
could well be of interest. 
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